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The objective of this study is to investigate the feasibility of a new pointing (position 
loop) controller for the NASA-JPL Deep Space Network antennas using the Disturbance 
Accommodating Control (DAC) theory A model that includes state dependent distur- 
bances was developed, and an example demonstrating the noise estimator is presented as 
an initial phase in the controller design The goal is to improve pointing accuracy by the 
removal of the systematic errors caused by the antenna misalignment as well as sensor 
noise and random wind and thermal disturbances Preliminary simulation results show 
that the DAC technique is successful in both cancelling the imposed errors and main- 
taining an optimal control policy 


I. Introduction 

Large, precision antennas for millimeter and submillimeter 
wave astronomical listening require precision pointing capabil- 
ity in the face of a host of nonlinear and random disturbances 
Included in this category of noise sources are structural mem- 
ber deflections under wind, thermal and gravitational loading, 
bearing friction torques, and hysteresis as well as electrical, 
optical, and mechanical misalignments introduced by sensors, 
thermal deformations, and structure model imperfections 
Traditional approaches for compensating systematic distur- 
bances rely on laboratory measurements and field data and 
employ open-loop (or feed-forward) compensation using static 
look-up tables to refine predicted target positions These tech- 
niques, although satisfactory in sub-X-band RF pomting, are 
marginal for the state-of-the-art telemetry requirements for the 
upcoming Voyager-Neptune flyby and beyond The augmen- 
tation of the deep space telemetry channel to provide Ka-band 
(32-GHz) capability to increase mission performance will 


require 1-mdeg pointing accuracy for feasible reception at 
distances greater than 20 AU The performance advantage 
between the current X-band and projected Ka-band is highly 
dependent on antenna pointing accuracy Successful deep 
space telecommunications will require the NASA-JPL 34-m 
and 70-m antenna pointing systems (see for example Fig 1) to 
exhibit pointing errors of 1 mdeg (rms) or better With current 
accuracy on the order of 5 to 10 mdeg, the antenna pointing 
loss at 32 GHz [1] as compared to loss at X-band (8 4 GHz) 
is magnified by the frequency squared The increased gam ad- 
vantage of the Ka-band could easily be lost without compar- 
able enhancement in the pointing accuracy, a performance 
requirement implicit to the higher gain antennas with narrower 
beamwidth 

Antenna pointing-tracking errors are typically functions 
of static and dynamic factors Mechanical misahgnment of 
sensors or inaccuracy in the predicts can be considered static 
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error sources, whereas dynamic factors would include wind, 
thermal and gravitational loadmg, etc The approach taken in 
this study is first to achieve a static error-free environment for 
precision pointing of the antenna The philosophy employed 
to achieve this goal is to treat the systematic misalignment 
errors, as well as the servocommands, as disturbances to the 
controlled system Hence, as an initial phase m the sequence of 
new controller design, this study addresses pointing-tracking 
in the presence of noise The research goal is intended to pro- 
duce the end product controller, hence, further investigation is 
necessary to augment the controller to include dynamic errors 
caused by random thermal and wmd disturbances 

Antenna pointing improvements can be developed through 
a sequence of progressive controller modifications solely in 
the existing software routines Current algorithms can be 
enhanced to simultaneously provide servotracking and correc- 
tion of the systematic errors as well as beam stabilization in 
the presence of random disturbance torques For a given 
antenna, servodnve, feed configuration, and surface distortion 
profile a computer software package could be developed to 
optimize the performance to achieve, adaptively, the maxi- 
mum antenna gain for a prescribed direction vector with a 
“smart” controller 

Typical antenna controllers consist of an analog rate 
loop and a position loop closed through a digital computer 
The control algorithms for the position control are either 
proportional-integral (PI) or state feedback control The PI 
control is accomplished by applying gains to the position 
error and the integral of position error The weighted sum of 
these signals is the commanded rate for the velocity loop 
Zero steady-state error to a ramp input is realized with the 
PI controller in this Type II system 

The more sophisticated method utilizing state feedback 
allows specification of the eigenvalues of the closed position 
loop The initial disadvantage of the state feedback is the 
requirement that all the states of the system be available 
for the feedback control The technique of state estimation 
has circumvented this problem, providing the controller 
with an estimate value for each of the unmeasurable or un- 
certain state signals The feedback gains are selected to yield 
the designer’s selected eigenvalues to achieve desired per- 
formance of the system This technique was incorporated 
in the upgrade of the 70-m antenna axis servos [2], [3] with 
the estimator gain vector selection based on system specifica- 
tions, minimal estimator error, and insensitivity to encoder 
and digital-to-analog (D/A) quantizations 

The unique controller enhancement proposed in this study 
suggests that, simultaneous with state estimation, another 


vector be estimated to represent the disturbance state which 
could then be used in determining a more complete control 
strategy The philosophy of the Disturbance Accommodating 
Controller (DAC) is based on the concept of state modeling 
of the disturbance vector as developed in [4] In the follow- 
ing section the basics of the DAC theory are discussed as 
applied to the antenna pointing problem A history of this 
DAC technique appears in [5] 

II. Theoretical Background 

The antenna pointing system can be modeled by a set of 
first-order linear differential equations 

x = Ax + Bn + Fw, x(0) = x Q ) 


where x(r) is an (n X 1) state vector, u(t ) is an (r X 1) control 
effort vector, y(t) is the (mX 1) output vector, and w(r) is a 
(p X 1) collective disturbance vector representing the com- 
bined effects of all the uncontrolled forces/torques acting on 
the antenna For instance, w(r) can be projected to include the 
effects of lateral winds, misalignments, etc It is assumed that 
w(r) cannot be predicted or measured accurately The matrices 
A, B, C, and F represent the continuous-time transmission 
matrices for the respective signals 

The state variable feedback is represented by the control 
law 

u = -K x x (2) 

where K x is a feedback gain matrix, not necessarily producing 
an optimum controller Two control cases will be analyzed fur- 
ther, the first assumes no disturbance F and the second case 
assumes that F exists 

A. Assuming No Disturbances, F = 0 

Substitution of Eq (2) into Eq (1) yields the closed-loop 
form with no disturbances, 

x = Ax - BK^x = (A - BK^)x (3) 

The advantage of the state feedback is the ease by which 
the closed-loop eigenvalues of the system, obtained from 
Eq (3), are arbitrarily specified through the selection of the 
gain matrix K x (also called the pole placement technique) 
However, the pole-placement method does not guarantee that 
the design is optimal On the other hand, if the optimal con- 
troller is designed, the quadratic performance technique, from 
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the theory of optimal control, provides the optimal steady- 
state solution to the minimum control effort and minimum 
transient deviation of the state from the origin problem, 1 e , 

u = -K g x (4) 

where the Kalman gam K y = -R - 1 B r P The term P is a sym- 
metric positive semidefinite solution to the steady-state 
Riccati equation in the matrix form, that is, 

A^ + PA+Q-PBR-'B^ = 0 (5) 


where e = x -x From Eq (8) it is apparent that the dynamic 
behavior of the error signal is determined by the observer 
poles If the matrix (A - Kq^C) is a stable matrix, the error 
vector converges to zero for any initial error e(0) 

Since A, C are fixed by the system, matrix Kq^ determines 
the estimator performance Again the pole-placement tech- 
nique can position the estimator eigenvalues from Eq (8) for 
proper performance, that is, x will converge to x regardless of 
the initial states x(0) and x'(0) Hence, the overall closed-loop 
noiseless system with full state estimator feedback can be ex- 
pressed in the state variable notation as 


The matrices R(X)) and Q(>0) are symmetric weighting 
matrices in the associated quadratic performance index for the 
continuous linear regulator, that is 


J(t 0 ) = x(T) t ?(T)x(T) 

and 


+ 



[x r Q(t)x + M 7 R(r)«] dt 


( 6 ) 




( 9 ) 


Note that the control is a time-varying state feedback, even if 
A, B, Q, and R are time-invariant, K g (t) varies with time 

Implementation of state feedback requires knowledge of 
the entire state vector In practice, however, not all state vari- 
ables are available for direct measurement Hence, a closed- 
loop estimator is utilized to predict the values of the unmea- 
surable state variables based on the measurements of the 
output and control variables The state estimator vector x' is 
described by 

x = Ax' + Bw + K q O' - Cx'), x'(0) = xL 


Note that the dynamics of the closed-loop system depend 
on the eigenvalues of both the controller and the estimator 
However, the separation principle allows the independent 
design of the controller and the estimator gam matrices assum- 
ing the observer poles are chosen correctly 

The optimal regulator described so far accommodates only 
initial conditions or impulse type disturbances and hence is 
incapable of tracking or handling typical noise inputs In the 
case of finite input disturbances, the control law of Eqs (2) 
through (9) cannot attain and maintain track, i e ,y(f) ^y c (t) 
where y c is the commanded output 


or 

x = (A-K 0 C)x'+Bn + K 0 y (7) 

u x u x 

where Kq^ is the estimator error gain matrix (n X m), and the 
eigenvalues of (A - Kq^C) are commonly called the observer 
poles Recall that the system must exhibit complete observa- 
bility in order to determine the state vector using the output 
and control variables 

The inaccuracy, e, in the state dynamics incurred in using 
the full-order ( n X 1) estimate x rather than the actual state x 
is given by subtracting Eq (7) from Eq (1) (with F= 0), l e , 

e = x - x = (A - K 0 C) e (8) 


B. Case of F # 0 

Consider now the plant equations in the form of Eq (1) 
with F ^ 0, 

x = Ax + Bu + Fw (10) 

where F is an ( n X p) matrix and w(r) is a p-dimensional dis- 
turbance vector Let us formulate the optimal regulator prob- 
lem in such a way that at the terminal time T, the resulting 
control law always brings the state x(t) and the velocity x(f) 
back to the commanded state and velocity, x c (f) and x c (t), 
respectively, in the presence of any finite constant disturbance 
w(t) = k With neither the noise nor the servocommand known 
a prion, treating w(f) as either a deterministic input or a non- 
deterministic input with a known probability is impractical 
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since reliable information about the disturbance is not avail- 
able The primary objective of the control is to manipulate 
u(t) in Eq (10) so that the output y{t) approaches and main- 
tains the commanded value y c (t) promptly Hence, the prob- 
lem is reduced to finding a control which minimizes the func- 
tional J(u), 


rr 

J(u) = lim / 
t-*Tj o 


[<* Qx ) + f(u, it, 


))dt (11) 


subject to the constraints imposed by Eq (10) where the 
vector disturbance w(0) = 0, and w(t) is assumed to satisfy 
the linear differential equation 


d p w 

~df 


+ a 


p - 1 


d p 1 w + 
dt p ~ l 


dw 

+a i~3r +a o w = 


0 


( 12 ) 


where the are known, real scalar constants 

The class of admissible disturbances w(t ) defined in Eq (12) 
can be characterized as the set of scalar functions 

w(t) = H z(t) (13) 

with H a real p X p matrix and where 

z(t ) = D z + a(t) (14) 

where z is a real p vector and D is real matrix (p X p) The a(t ) 
in Eq (14) represents the uncertainty in the noise model This 
representation of w(r), illustrated in Fig 2, shows that the 
optimal controller is designed by first building a duplicate 
model for the disturbance process typified by Eqs (13) and 
(14) This noise estimator is driven by the vector Cx(t), as is 
the plant state estimator The noise estimate and the system 
state estimate are weighted and summed to yield the control 
law u(t ) As t-*T, x (t) approaches the steady -state x(r) =x c (t) 
prior to any change in the state command 

Mathematically the noise is not precisely known Hence, 
the a(t ) represent completely unknown sequences of random- 
intensity, random-occurring, isolated delta functions Antenna 
experimental data have shown that the alignment uncertainty 
exhibits the less than noisy properties of a stochastic process, 
and thus, in this particular case, Eq (13) appears to be a rea- 
sonable model of the systematic disturbances 

The control law is effectively divided into two parts, that is. 


u „ + u. 
s a 


( 15 ) 


where u s is assigned the task of servoregulation, and the com- 
ponent u d is responsible for counteracting both disturbances, 
w(r) and the reference signal y c (t) m servotracking In the case 
of regulation, y c (t ) = 0, the control u d must be capable of 
satisfying the relation 

B(f)« d (0 = -F(f)w(0 = -FHz(r) (16) 

for all admissible w(t) in order to cancel out the disturbance 
F Johnson in [5] has shown that if the rank of [B|FH] = 
rank of [B] , then 


-F(0H(0 = -B(f)K z (r) (17) 

for some gain matrix K z (f) Failure of this condition means 
that some residual effect of w(f) will always occur Assuming 
Eq (17) is satisfied, the counteraction torque, u d , is 

u d (t) = -K Z {t)z{t) (18) 

and the regulation control, u s , is chosen by the conventional 
means as 


«,(0 = (19) 


The open-loop system augmented to include the noise state 
is described by 
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( 20 ) 


The exact closed-loop state and disturbance state vectors, 
assuming all states are available for measurement, are described 
by the differential equations 


X 


A - BK, 

FH - BK, 

z 


0 

D 



( 21 ) 


Redefining the representation of Eq (20) so that the new 
state vector includes both system and disturbance states gives 
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Ax+Eu+8 


where K 0 , the composite estimator gain matrix, is chosen to 
force x'(f) -* x(t) Using the same principle mentioned earlier 
in Eq (8), the eigenvalues of (A - KqC) are chosen for the 
augmented system The practical control law for servoregula- 
tion in the face of noise becomes 


y = Cx 

Examination of Eq (21) reaffirms the criterion set forth in 
Eq (17) to eliminate the noise term in the system state equa- 
tions However, an obvious problem arises how is it possible 
to access the actual disturbance state z(r) 9 Of course, z(t) is 
not completely measurable, but it is possible to resort to esti- 
mator theory once again to observe and predict the noise state 
in an approach similar to that used in the plant state estima- 
tion Hence, an appropriate control is obtained by replacing 
the actual noise state in Eq (18) with the estimate of the noise 
state le , where estimates of z(f) and x(t) can be ob- 
tained from y(t) by on-line, real-time state reconstruction 
In general K z (r) is shown to be not unique [4] 

The actual closed-loop plant state x(r), with the assumption 
that all states are available for measurement, and the disturb- 
ance state error variable e z (t) are described by the differential 
equations 


i: 


~A-BK x 

FH - BK 0z " 


X 

+ 

0 " 

Jz_ 


0 

D + FH 


e z 

L J 

o(t) 


nit) = -Kx(t)x'(t) = -K^x'^-K (r)z'M (25) 

where z', x are the estimates of z(t) and x(r), respectively, 
obtained from y(t) by on-line, real-time state reconstruction, 
and -K^OxV) ls the control required to minimize a perfor- 
mance index if the disturbances were not present m Eq (10) 
In general K z (r) is not unique, as shown in [4] 

The system is reorganized now to include servotracking 
rather than just servoregulation Similarly the servocommand 
can be treated as a “disturbance” to the plant Recall that the 
primary control objective in the antenna pointing system is 
that of servotracking .y c (f), where in general the command is 
related to the system variables (xj , , x n ) by the equation 

y c it) = C(t)x(t) (26) 

In tlmcase, the objective is to control the plant output^ (f) so 
that C is equal to C in Eq (1) The behavior of y c is assumed 
expressible by the servocommand model 

y c it) = G(r)c(r) 

(23) c = E(r)c(r) + Ai(r) 



where e = z - z', with Kq z the appropriate observer gain for 
the disturbance state estimator 

Examination of Eq (23) reveals that if e z is zero, the be- 
havior of x(r) is totally independent of the disturbance w(t) 
Of course, the ideal case is rarely realized, however, Ko z (t) 
is chosen so that e z (/)-K) rapidly for all initial values of x 0 , 
z 0 , e 0 so that the closed-loop plant state x(i) is essentially 
insensitive to external disturbances that can be generated by 
Eq (14) 

Using the composite model of Eq (22), the corresponding 
augmented state estimator vector x'{t) is described through 
observer theory as 


where G(r), E(r) are determined beforehand by appropriate 
modeling procedures, and c represents the servocommand 
state vector The vector j a(t) represents the uncertain impulse 
sequences, similarly introduced in the disturbance model in 
the form of o(t) Note that in the case of set point regulation, 
y c is essentially a constant, and hence E(r) = 0 and G(r) is 
the identity matrix with the assumption that the y t are inde- 
pendent outputs In servotracking, the y c (t) are allowed to 
vary continuously with time and E(f) is chosen accordingly 
Exact servotracking cannot be realized unless the servocom- 
mand error e = Gc - Cx is zero Hence, u s must be chosen so 
that (y c -+y) rapidly approaches the null space of C for all 
initial conditions 


or 


x' = A x' + B u + K 0 (2 - C x') 


x = (A-KqCJx' + Bu + K^ 


(24) 


Hence, the three individual plant, disturbance, and servo- 
command models can be combined into a single composite 
open-loop model 


Ax + B u + 6 


(28) 
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where 



A 

FH 

0 


B 


0 

A = 

0 

D 

0 

, B = 

0 

, i = 

o(t) 


0 

0 

E 


0 


tit) 


y = C x where C = [ C | 0 | 0 ] 


Figure 2 illustrates the continuous-time system with the 
DAC installed The solution approach is as follows An appro- 
priate generic control law for servotracking m the face of noise 
follows as 

u = -N(r)z'(0 - K^Ox'Cf) - K c (t)c'(0 (29) 

Utilizing the servocommand state estimator, the closed-loop 
plant can be described as 

x = (A - BI^)* + (-BK 2 + FH)z' - BK/ 

+ B [ -K z (z - z') - K x (x - x') - K (c - c') ] (30) 

With the appropriate choices of the gams, K x , K z , K c , the esti- 
mate errors will approach zero quickly and the noise terms in 
z should have no effect on.y(0 


ous model, is required for the investigation Reference [6] 
derives the discrete composite state variable representation 
with the assumption that the noise is not state dependent, i e , 


*[(ii+1)7’]' 


'A(nT) 

FH(nT ) 
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~ x (nT) 

z[(«+ IF] 
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D(nT) 
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z(nT) 

_c[(«+i)n 
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E(nT) 


c(nT) 



B(nT) 


y(nT) 

+ 

0 

u(nT) + 

o(nT) 


0 




x(nT) 


y(nT) = [C(«r)|0|0] 


z(nT) 


c(nT) 


(31) 


where the discrete-time plant matrices are 


III. Mathematical Models for the 
Discrete-Time DAC 

The continuous-time model of the disturbed dynamical sys- 
tem described m Eq (1) can be transformed into a discrete- 
time prototype for investigation of an analog plant controlled 
by a digital computer For simplicity, assume that the signal 
sampling time is coincident with the control effort application 
time In feedback controller designs, control decisions are de- 
termined in light of real-time data provided to the controller 
through sensors The term “sampled-data” denotes data pro- 
vided to the controller which are updated only at specific 
isolated points in time Between these updates, the data pro- 
vided to the controller typically are held constant Likewise, 
the control decisions are updated only at specific isolated 
times In between the decision updates, the control action 
u(t) either remains constant or follows a prescribed interpola- 
tion rule The computer or digital controller is capable of pro- 
cessing only sampled-data and executing discrete-time control 
policies usually written as difference equations 


An appropriate discrete-time representation of the system, 
disturbances, and servocommands, analogous to the continu- 


A(nT) = <b A [[f 0 + («+l)7\r 0 +«r]] 


B(nT) 


4 plant state transition matrix in discrete time 
H)T 


7[(«+on = 


/■f 0 +(«+i 

J, 0+ nT 
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$4 [t 0 + (n + 1)7 \t] B(t) (It 


r 0 +(« + l)r 


^[r 0 + (n+l)T,r] F(r)H(r) 


Jt+nT 


dr 


,(TX)o(m 

+n T 
rt 0 +(n + l)T 

FH(nT)= / * A |[f 0 + («+l)7>] 

J t o +nT 

XF(t)H (T)* D (T,t 0 + nT)dT 
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and the discrete-time noise terms become 


where 


D{nT) = % [[r 0 + (#i+l) 77 0 + «rj] 
4 noise state transition matrix 


a(nT) 


J ^r 0 +(n + i)r 
l o +nT 


<h D [t 0 +(n+l)7-,?j o$)<$ 


(32) 


In the time-invariant case, the matrices A, B, C, F are con- 
stant element matrices, and Eq (32) is simplified to 


E{nT) = % [(« + \)T/iT] 


4 discrete transition matrix for 
the servocommand 


u(nT) 


fV" + 

Jt n +rxT 


1 )T 


X<t> E [t 0 + (n+\)T,$] 


( 35 ) 


A = e AT 


D = e° T 


B= e A(r-T) Bdr 
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e A (7'-r) FH 


I. 

/; 

X ^ e D ( T ~U ott + t 0 +nT)d$ 


FH = / e MT - T) FH e Dr dr 


I. 

-s: 


(P<- T -n o(t + r +fi7 , )rff 



Similarly, the discrete-time servotracking state 
Eq (27) can be represented in the form 

y c (nT) = G{nT)c(nT ) 

£•[(« + i)r] = £(nr)c(«r)+j«(/ir) 


model of 



In the time-invariant case, Eq (35) is reduced to 


E(nT) = e ET 


1 , 


d(nT)= / ju« + t 0 + nf) 


(36) 


The model may be generalized further to include various 
exceptional case studies [5] For example, the antenna sys- 
tematic errors appear to be dependent on the particular 
azimuth/elevation position of the target, hence, the noise 
w(t) can be made a function of the system orientation or 
of the state of the plant In this case the disturbance model 
can be augmented to include the state dependency by adding 
extra terms as follows 


w(nT ) = H(nT)z(nT) + L(nT)x(nT) 


z[(n + 1)7’] = D(nT)z(nT) +M(nT)x(nT) 
+ a(nT ) 


(37) 


In this case it is necessary to derive the appropriate rela- 
tionships from the continuous to discrete-time case in order 
to ascertain the mathematical meaning of the additional terms 
L(nT ) and M(nT) in Eq. (37) These relations are given in 
Eq (38) 
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(39) 


where in the time-variant case, 


where u(nT), y(nT) and y c (nT ) denote the inputs to the ob- 
server The matrices Kq x , K 0z , Kq c are arbitrary gain matrices 
that the DAC designer selects in accordance with the desired 
estimator response 


7 2 K« + i)n 



e MT-T) FH 


X 



e D(r_t) x(£ + t 0 +nr)d£ 


dr 


In order to establish the dynamics of the estimators, con- 
sider the state errors e x , e z and e c defined as 

e x = x(nT)-x'(nT) \ 

e z = z(nT)-z’{nT) > (40) 

e c = c(nT)-c'(nT) ] 


and 


fT 

M(nT) = / e D(r - f) L x(t + t Q + nT)dr 

Jo 


Construction of the on-line, real-time estimation of the 
three states-plant state x(t), disturbance state z(r), and the 
command state c(f)— requires the discrete-time state estimator 
given by 


These error dynamics are described by the discrete equations 
derived from Eq (39), considering only first -order variations, as 
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In order to produce reliable estimates, the observer gain 
matrices (K 0x , K Qz , K Qc ) are chosen so that the errors in 
Eq (40) decay toward zero rapidly between control updates 
Hence the homogenous solution of Eq (41) is made asymptot- 
ically stable to the errors equal to zero In general, Eq (41) is 
a time-varying set of difference equations The gains K 0x , 
K 0z , K 0c can be solved using the discrete Riccati equation 
from optimal control theory For the case of constant element 
matrices A, C, FH,D , the design of the estimator gain matrices 
can be accomplished by the conventional eigenvalue placement 
method Defining the error dynamics as a system with a char- 
acteristic matrix A e 

'(A+K 0x ) FH 0 

A g = K q C D 0 (42) 

0 0 ( E + K 0 G ) 

the eigenvalues are positioned suitably (say, at zero) within 
the unit circle A block diagram of the composite observer is 
shown in Fig 4 


IV. Determination of the Discrete-Time 
Control Function u 

The determination of the control function u in the discrete- 
time case involves several steps 

(1) The state estimators are weighted and summed to 
determine the control law, l e , 

u{nT) = f[x'(nT), c'{nT),z'(nT), nT] (43) 

(2) The control function is divided into two subtasks as 
mentioned previously, 

u(nT) = u s (nT) + u <] (nT) (44) 

where the component is responsible for the servoreg- 
ulation and the u d effort is assigned the task of distur- 
bance removal including servotracking Substitution of 
Eq (44) into Eq (27) yields the plant state relation 

*((«+ l)r) = A(nT)x(nT) + B(nT)u s (nT) 

+ B{nT)u d (nT) + FH(nT)z(nT ) 

+ y(nT) (45) 


Since the control effort in a discrete-time control problem 
is usually held constant between two consecutive sampling 
times, it is impossible generally to remove all the disturbance 
effects Likewise, the presence of the uncertainty sequence 
y(nT) also limits the idea of complete time cancellation of the 
noise Hence, the concept of “complete cancellation” means 
only that the noise effects FHz(nT ) are removed as they 
appear at isolated sample times, i e , 

B(nT)u d (nT) + FH(nT ) z(nT) + E(nT) c(nT) =0 

(46) 

The condition for existence of u d (nT ) to satisfy Eq (46) is 
[FH I A e ) = B[K z | K c ] 

Complete disturbance cancellation exists if, and only if, 

FH(nT ) = -B(nT)K z (nT) 
and 

E{nT) = - B(nT)K c (nT ) (47) 

for some matrix K z {nT) and K c (nT) Assuming the conditions 
of Eq (47) are satisfied, the control u d {nT) can be chosen in a 
practical sense as 

u d (nT) = ~K z {nT)z\nT)~K c {nT)c\nT) (48) 

where z\nT) is the noise state determined by on-line, real- 
time estimation of z(nT ) The enclosed loop error dynamics 
using e z , e c may be incorporated into the model as 

x [{n + 1)7"] = A{nT) x(nT ) + B(nT ) uJnT) 

- [B(nT)K z (nT)} e z (nT ) 

- [B{nT)K c (nT)} e c (nT) + y{nT) 

(49) 

Hence, the noise effects have been reduced to the ( BK z e z 
+ BK c e c ) term, which should decay rapidly toward zero, and, 
of course, the isolated uncertainty sequence y(nT) The servo- 
regulating control, u s (nT), can now be designed by conven- 
tional methods assuming the noise has been removed 

A complete block diagram of the original continuous-time 
plant model, and the proposed DAC with full state discrete- 
time composite observer, is shown in Fig 4, with the control 
law 

= -xj 
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and 


(56) 


u = -K x 

S X 


u d = -Kz -Kc (50) 

V. Summary of DAC Procedure 

(1) The disturbance w(t ) is determined experimentally to 
ascertain distinguishing characteristics Suppose that 
w(f) is noted to consist of an uncertain bias at times 
and in other intervals w(r) exhibits uncertain ramp fea- 
tures Thus, the waveform of the disturbance has a 
general form 


The control law u d has the form 

= ~K z ~ K c c ( 57 ) 

where K z and K c are chosen to ensure disturbance re- 
moval and tracking, respectively 

(4) With available software tools, the equations can be 
incorporated into a single system simulation-see 
Eq (28)-and tested for various noise inputs and 
servocommands 


w(t) = k 1 +k 2 t (51) 

where k l , k 2 are unknown constants which change 
value at unknown times 

(2) With the description of w(r), the designer determines 
the simplest differential equation model for this class 
of disturbances, that is, the lowest order differential 
equation for which Eq (51) is the general solution 
The corresponding DAC matrices H, D are determined 
from the general form given in Eqs (13) and (14) In 
observable canonical form the model becomes 


VI. Test Model 

A simplified test model for an antenna servomechanism is 
used to illustrate the DAC procedure (refer to Fig 3) The 
objective is to control the elevation of an antenna designed 
to track an RF signal The antenna and its drive mechanism 
have a moment of inertia J and damping B arising from bear- 
ing friction, aerodynamic friction, and the back emf of the 
dc-servodnve motor The equations of motion are 

J9+B9 = T c + T d (58) 


w(t) = [ 1 



where T c is the net torque developed by the drive motor, and 
T d represents the disturbance torques possibly due to wind, 
static misalignments, etc Substitution of the assumed coeffi- 
cients in Eq (58) yields 
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where (& b 2 ) are constants The characteristic poly- 

nomial of the disturbance model D is equated with the 
resulting characteristic polynomial attained from the 
Laplace transform of Eq (51), that is, 

det[Xl-D] = \ 2 +b 2 X + b 1 = X 2 (54) 

and in this example, b x and b 2 equate to zero 

(3) The controller problem is separated into two subtasks, 

u = u s + u d (55) 


£+10 + 6 9 - u + w (59) 

In this example, the coefficients for Eqs (58) and (59) 
were arbitrarily selected and may be unrealistic They were 
selected, however, to describe the effect of the new DAC con- 
troller The general shape of the servocommand angle 6 c {t) is 
assumed to be composed of step and ramp functions Hence, 
the servocommand 9 c (t) is modeled by y c estimated using the 
commanded rate 0 c (f), and the acceleration 9 c (t) For the pur- 
pose of maintaining a good tracking accuracy, it is reasonable 
to assume that the antenna drives are capable of following the 
peak velocity 8 c (t ) in the steady state with acceptable error 
Since the objective is to permit acceptable communication 
signal reception, the dependence of the signal amplitude on 
pointing error is a major concern The corresponding servo- 
command can be represented by the following state-space 
representation 


such that u s performs servoregulation and u d is respon- 
sible for servotracking a command input and the distur- 
bance removal Hence, the control effort u s is chosen 
for regulation assuming no noise, typically, 


c(f) = E(f) c(t) + n(t) 
y c (0 = G(r)c(r) 


(60) 
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Similarly the disturbance w(t) (assumed to be step and ramp 
torques in this example) can be modeled as suggested in 
Eq (54), 


w(t) = H z(r) 
z(t) = D z(t) + o(t) 


(61) 


Thus, the open-loop system can be represented in the form of 
Eq (28) where 



'o 

f 


'0' 

A = 

-6 

- 1 _ 

F = B = 

_ 1 _ 


'o 

r 



E = 

0 

0 

G = C = [ 

1 


H = [ 1 0 ] 


The state feedback gains have been chosen with K x = [94, 
19] to ensure the undisturbed system closed-loop poles at 
-10, -10 The estimator poles for the plant, noise, and servo- 
command states are selected in the usual manner using pole- 
placement techniques for approximately 3 to 5 times faster 
response than the combined plant and feedback controller 
A computer simulation of the closed-loop model in Eq (28), 
shown in Fig 5, is used to demonstrate the effectiveness of 
the DAC in disturbance rejection Figure 6(a) illustrates the 
controlled output variable 0(f) servotracking the command 
input 0 c {t) without the DAC The disturbance assumed in 
this example is plotted in Fig 6(b) with the controlled vari- 
able y(t) Without control other than state feedback, the 
output is unable to distinguish the control input from the 
disturbance and tends to follow the noise signal rather than 
the servocommand 


In Fig 7(a) and (b), the same example with the inclusion of 
the DAC demonstrates the effective servotracking of the com- 
mand m the presence of the noise input Note that only a 
slight perturbation occurs in the controlled variable y (f) at 
approximately 5 0 sjust as the disturbance has occurred 

The analogous discrete-time system has been simulated to 
demonstrate the degradation expected in tracking when the 
position-loop of the controller is implemented via a digital 
computer Figures 8 and 9 display the discrete-time system 
in a noisy environment both with and without the DAC incor- 
porated in the loop, assuming a sampling time T- 0 1 second 


VII. Conclusions 

The feasibility of implementing a disturbance accommodat- 
ing controller has been investigated as applied to an analog 
servodnve for positioning an RF antenna The DAC is designed 
for synthesizing and rejecting waveform-structured distur- 
bances The form of the systematic pointing errors inherent 
in antenna tracking systems appears viable to this character- 
ization of the disturbance as structured waveforms rather 
than the noise generated through random processes with sta- 
tistical descriptors The waveform type of disturbances can 
be modeled according to a prion data by determination of the 
corresponding differential equation, and hence, the state 
representation of the waveform structured noise 

In this study, simulation results show that the DAC is an 
appropriate technique for cancellation of the systematic 
errors, while simultaneously allowing an optimal control 
policy to regulate the system The ease with which the DAC 
is implemented along with the existing servo-control is another 
attribute of this technique Practical implementation issues 
such as model order, computation time, and storage require- 
ments offer no expected challenges for microprocessor-based 
controllers Further study is necessary to incorporate the 
state-dependency issue in regard to systematic pointing errors 
expected in antenna position controllers 
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Fig 1 Overview of the 70-m antenna pointing system at JPL 
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Fig 5 Simulation model for the DAC servocommand example 
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OUTPUT 



Fig 6 Servocommand response for the continuous-time example 
without the DAC Installed, (a) servocommand and controlled 
variable; (b) disturbance input and the controlled variable y(t) 



Fig 7 Servocommand response for the continuous-time example 
with the inclusion of the DAC. (a) servocommand and controlled 
variable; (b) disturbance input and the controlled variable y(t) 
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Fig. 9. Servocommand response for the discrete-time example 
with the DAC installed, (a) servocommand input and controlled 
variable; (b) disturbance input and the controlled variable y{t) 


Fig. 8 Servocommand response for the discrete-time example 
without the DAC installed (a) servocommand input and controlled 
variable, (b) disturbance input and the controlled variable y(t) 
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